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ABSTRACT 

\^ • We investigate the potential of exploiting the Sunyaev-Zeldovich effect (SZE) to study 

(— I ' the properties of hot gas in galaxy groups. It is shown that, with upcoming SZE sur- 

veys, one can stack SZE maps around galaxy groups of similar halo masses selected 
I ' from large galaxy redshift surveys to study the hot gas in halos represented by galaxy 

O I groups. We use various models for the hot halo gas to study how the expected SZE sig- 

nals are affected by gas fraction, equation of state, halo concentration, and cosmology. 
Comparing the model predictions with the sensitivities expected from the SPT, ACT 
and Planck surveys shows that a SPT-like survey can provide stringent constraints on 
the hot gas properties for halos with masses M >, 10"'^'' h~^MQ. We also explore the 
' idea of using the cross correlation between hot gas and galaxies of different luminosity 

to probe the hot gas in dark matter halos without identifying galaxy groups to rep- 
resent dark halos. Our results show that, with a galaxy survey as large as the Sloan 
Digital Sky Survey and with the help of the conditional luminosity function (CLE) 
model, one can obtain stringent constraints on the hot gas properties in halos with 
masses down to 10^"^ h~^MQ. Thus, the upcoming SZE surveys should provide a very 
promising avenue to probe the hot gas in relatively low-mass halos where the majority 
of L* -galaxies reside. 
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1 INTRODUCTION from the hot phase, various feedback processes, such as su- 

pernova explosions and AGN activity, may act as an efficient 

In the current cold dark matter (CDM) paradiem, the laree- r 1 .■ . 1 ■ .1 ■• r,i 

^ y r o ) o source 01 reheatmg, strongly impactmg the properties ot the 

scale structure m the Universe grows due to gravitational , . „ ^1 ^ r 1 ^ 1 1 1 

, r 1 , 'lot gas. Consequently, the amount ot hot halo gas and its 

instability, forming dark matter halos withm which galaxies 1- , .1 ,• 1 i-rr- . r , j ^ r 

, distribution may be very dmerent trom that expected trom 

form. Since the cosmic gas is initially cold and well mixed . , ,. 1 , t 1 1 .1 , , 1 . r 

the simple non-radiative model. Indeed, the total amount ot 

with the dark matter, the gas component is expected to fol- , . . , , , , 1 i- ■ 1 1 

, hot gas m low-mass halos, such as those hosting spiral galax- 
low the evolution oi the dark matter component until it is 

heated by accretion shocks during halo formation, resulting 



ies, is observed to be much lower than that expected from 

the universal baryon fraction. Even for clusters of galaxies, 
m extended gaseous halos that trace out the potential wells , j_i 1 ^ r. ^- ■ ^ w 1 1 1 ^ 

° . where the hot gas fraction is observed to be comparable to 



of dark matter halos. In the simplest case where the pro- 
cess is assumed to be non-radiative, the total amount of hot 
gas contained in a gaseous halo is expected to be roughly the 
amount of gas initially associated with the dark matter halo, 
and the distribution of the gas is governed by hydrostatic 



the expected universal fraction, the state and distribution of 

the hot gas appear to be strongly affected by star formation 

and AGN activity. Thus, the current state of hot gas in dark 

matter halos reflects not only the outcome of gravitational 

. ^ ,. , , . . . , collapse, but also contains important information regarding 

equilibrium. In reality, however, the situation IS much more , „ ,. , , ,. , 1 , 1 

galaxy formation and evolution withm dark matter halos. 
complicated. In addition to radiative processes which cause „, „ , , , , r . 1 1 . 1 1 ■ 1 ■ 

rr , therefore, the study ot the hot halo gas can provide impor- 

the halo gas to cool and form stars, effectively removing it , , , . , , r , ■ 1 1 . • • 

tant constraints on galaxy formation and evolution m a way 

complementary to that provided by observation of stars and 
* E-mail:Iir@bac. pku.edu.cn cold gas. 
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The hot halo gas has traditionaUy been observed 
through its diffuse X-ray emission. However, since the 
X-ray emissivity depends sensitively on the gas density, 
this method is only powerful for hot gas with relatively 
high density, such as that in the inner parts of massive 
clusters/groups of galaxies. An alternative way to probe 
the hot halo gas is through the Sunyaev-Zel'dovich effect 
(hereafter SZE) the hot gas generates in the cosmic mi- 
crowav e background (CMB) due t o inverse-Compton scat- 
tering (|Sunyaev fc ZeldovichlHoT^ '). Because CMB photons 
gain energy as they traverse hot gas, the surface brightness 
of the CMB (in a given band) in the direction of a clus- 
ter/group will be changed. In the Rayleigh- Jeans tail (where 
hpv/kB <^ TcMB ~ 2.73K, with fce the Boltzmann constant 
and hp the Planck constant), the change is a decrement. 
This observable decrement is proportional to the electron 
pressure intergrated along the line of sight, and is more sen- 
sitive than X-ray emission in probing hot gas within regions 
of relatively low density, such as in the outskirts of dark mat- 
ter halos. Moreover, the SZE "brightness" (i.e. the change in 
CMB brightness due to the SZE) is independent of redshift, 
making it possible to study the evolution of hot gas in dark 
matter halos out to high redshift. 

The first detection of the thermal SZE with high 
signifi cance was reported in 1978 (|Birkinshaw et ahl 
19781'). six years af t er th e concept was proposed by 



Sunvaev fc ZeldovichI (|l972l '). Since then, technological ad- 
vances have been such that this method has become an 
importan t tool in probing the hot gas in g alaxy clus- 
ters (e.g. iMvers et all Il997l : ICrego et all I2001I : iNord et all 
I2OO9D . As mentioned above, since the SZE surface bright- 
ness is independent of redshift, it is also considered a pow- 
erful tool to detect and study galaxy clusters at relatively 
high redshift, which can be used to con strain the rate of 
structure formation in the Universe (e.g. Birkinshawl 1 19991 : 
iHolder fc CarlstromI [2OO1I : iMuchovei et al.l |2007| ). In addi- 
tion, the SZE of galax y clusters ha s been combined with 
X-ray observation (e.& Turned 20021: LaRogue et al.l 2006 ) 
and weak lensing (e .g. lBartelmannll200ll : ISealfon et al. 
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lUmetsu et al.ll2009l ) to accurately probe the shape and mass 
distribution of clusters of galaxies. 

Despite this progress, the application of the SZE is still 
in its early stage, and the observational samples available to- 
day are still small. However, this situation will change dras- 
tically in the near future. New ground-based telescopes such 
as the South Pole Telescope (SPTj3 and the Atacama Cos- 
mology Telescopy (ACTfl as well as the Planclfl satellite, 
are expe cted to detect thousa nds of clusters through the 
SZE(e.g. iMoodlev et al.l [2OO8I ). However, even with these 
new observations, only the hot gas profile of the rela- 
tively massive systems (with masses larger than a few times 
10^"* /i~^M0) are expected to be studied individually. In 
the case of less massive halos, which host the majority of 
bright galaxies in the Universe, the data will be insuffi- 
cient for individual detections. However, one may bypass 
this restriction by stacking large numbers of groups to- 
gether in order to increase the signal-to-noise (hereafter 



^ http://astro.uchicago.edu/spt/ 
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S/N). Such an analysis, however, requires a pre-selected 
sample of groups. Fortunately, large surveys of galaxies 
such as the 2-degree Fie ld Galaxy Redshift Survey (2dF- 
GRS; ICoUess et ahl I2OOII) an d the Sloan Digital Sky Sur- 
vey (SDSS; York et al.l 2000| ) are now available for select- 



ing large and uniform samples of g alaxy groups that cover a 
wide range of halo masses (e.g.lCoto 2005; Mille r et al.ll2005 



iBer lind et all I2OO6I : lYang et all [iooa 120071 : IKoester et all 
I2007i r This makes it possible to study the average prop- 
erties of hot halo gas in relatively low-mass halos. 

In this paper, we explore the potential of using stacks 
of galaxy groups to probe and study the hot gas in relatively 
low-mass halos through its SZE. We base our investigation 
on groups selected from current galaxy redshift surveys com- 
bined with the expected SZE data from on-going surveys 
such as Planck, ACT and SPT. We consider a number of 
plausible models to describe the properties of the hot halo 
gas, and compare our predictions with the detection limits 
expected from ongoing surveys. In addition, we also explore 
the possibility of probing the hot gas in low mass dark mat- 
ter halos by using the cross-correlation between SZE maps 
and galaxies of different luminosities, which has the advan- 
tage that it does not require a pre-identification of galaxy 
groups. 

The structure of this paper is as follows. In Section [2] 
we present our models for dark matter halos and for the hot 
halo gas. In Section [3] we describe the SZE expected from 
dark matter halos of given masses. Section |4] describes in 
detail how to calculate both the cross-correlation between 
hot gas and dark matter halos, and between hot gas and 
galaxies of a given luminosity. In Section (5] we discuss the 
properties of the group catalog to be used in our modeling. 
Section |5] contains our predictions of various models and 
their detectabilities with the on-going SZE surveys such as 
SPT, ACT, and Planck. We summarize our conclusions in 
SecE] 

Unless specified otherwise, we adopt a ACDM cos- 
molog y with parameter s given by the WMAP 3-year 
data (Spergel et al.ll2007l . hereafter 'WMAP3 cosmology'): 
^^,0 = 0.238, Oa.o = 0.762, h = 77o/(100kms-^ Mpc'^) = 
0.734, and aa = 0.744. 



2 DARK MATTER HALOS AND HOT GAS 
2.1 Dark matter halos 

We assume that d ark matter halos follow the NFW 
l|Navarro et al.ll 19971 ) profile. 



^ a X ^ c 

else 



Here x = r/r^, with Vs a scale radius related to the halo virial 
radius r^ir via the concentration parameter, c = rvir/f's, and 
Pcrit is the critical density of the Universe. The characteristic 
over-density. So, is related to the critical over-density of a 
virialized halo, Avir, by 



So = 



3 ln(l -I- c)-c/(H-c) 



(1) 



http: / / astro, estec. esa. nl /Planck / 



For the ACDM cos mology considered here, we adopt the 
form of Avir given bv lBrvan fc NormanI l|l998l ) based on the 
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spherical collapse model: 



IStt^ + 82[nm(z) - 1] - 39[nn,{z) - 1] 



(2) 



where Q,in{z) is the cosmological density parameter of mat- 
ter at redshift z. Thus, the concentration parameter c is 
the only parameter required to specify the density pro- 
file of a halo of mass M. Numerical simulations show that 
halo concentration decreases gr adually with halo mass (e.g. 
iBuUock et all I2OOII : lEke et al.l [2OOL '). However the exact 
mass-dependence of the concentration parameter has not 
yet been well constrained, neither theoretically nor observa- 
tionally. Various fitting formulae based on nu merical simula- 
tions have been proposed in the literature (e.g. Bullock et al.l 



200ll:IZhao et al.ll2003l : iDolag et al.l2004l : iMaccio et al.ll2007l : 



Zhao et al.ll2009l '). For the purpose of this paper, the differ- 
ence between these different formulae does not affect our 
results qualitatively. In what foll ows, we adopt the fi tting 
formula for c = c(M, z) given by iMaccio et al.l (|2007l '). For 
simplicity we do not include the scatter in c for halos of a 
given mass at a given redshift. 

2.2 The Hot Gas Component 

By far the most common model used to describe the den- 
sity distribution, Pg{r), of hot gas in clusters is the /3- 
model, which consists of a constant density core and falls 
off as pg oc at large radii. Recent observation, however, 
has shown that the /3-model is not a good approximation 
for the gas density distribu tion outside the core regions of 
clusters (e.g. Vikhlinin et al. 2006; Komatsu & Scljak 2001; 
iHallman et all 120071 ). In this paper, we adopt a model in 
which the gas is assumed to be in hydrostatic equilib- 
rium(HE) with a given equation of state. This assumption is 
supported by numerical simulations (c.g.' Evrard et al]|l996l : 
iBrvan fc Norman 1998: Thomas et al. 200ll), and has as an 
additional advantage that it is straightforward to link the 
gas properties to those of the dark matter halos. 

For an isothermal equation of state, the virial temper- 
ature of a halo is defined as 



pm^ GM 
2 fce rvir 



(3) 



where nip is the proton mass. We assume a fully ionized gas 
with a mean molecular weight p = 0.588. In reality, however, 
the hot gas within dark matter halos is not expected to be 
isothermal. Therefore, we consider two alternative models; 
one in which we assume a more realistic, polytropic equation 
of state, and one in which we consider some amount of initial 
entropy injection. In either case, the density profile of the 
gas is obtained by solving the hydrostatic equation. 



d^ 
dr 



Gpg(r)M(< r) 



(4) 



where M{< r) is the total mass within radius r. For sim- 
plicity, throughout this paper we will ignore the gravity of 
the gas. 

2.2.1 Polytropic Model 

Let us first consider the polytropic case, which we regard as 
our fiducial model. The equation of state in this case is 

Pg oc , (5) 



with r the polytropic index. Throughout this paper we 
adopt r = 1.2, in agre ement with both observations (e.g. 
iFinoguenov et al. 2007) and with numerical simulations 
fe.g. iLewis et al. 2000: Borgani et al.|[2004l : lAscasibar et al.l 



With this equation of state, it is straightforward to 
obtain the density profile pg and the temperature profile Tg: 



and 



Psi^) = Pg,o??poiy(a;) , 



(6) 



(7) 



where, as before, x = r/vs, and Tg^o and pg,o are the central 
temperature and central density, respectively. The function 
?7poiy is dimensionless and is obtained from the hydrostatic 
equation: 



T?poly(3::) = [1 -~Bg{x)] 



i/(r-i) 



where 



and 



with 



g{x) = 1 



ln(l + x) 



(8) 



(9) 



B = 2 



c 

m 



fix)=ln{l + x) 



1 + x 



(11) 



For given halo mass and redshift the density and tempera- 
ture profiles of the gas are therefore specified by three pa- 
rameters: r, Tgfi and Pg,o- We set F — 1.2, and the other 
two free parameters are specified by boundary conditions. 

As the first boundary condition, we assume that the 
temperature at the virial radius is equal to the virial temper- 
ature, in good agreement with num erical simulations (e.g. 
iFrenk et al]|l999l : iRasia et all 120041 '). The central tempera- 
ture of the polytropic gas can then be written as 



1 + 2- 



1 cg{c) 



/(c) 



(12) 



To specify Pg,o, we fix the total gas inside the halo to be a 
universal fraction of the halo mass. In particular, we assume 
that a fraction /star of the gas initially associated with the 
halo either has formed stars or is in the cold phase. The 
central density of the hot gas is then given by 



Pg,0 = (1 - ./star) 



fib 



fin 



47rr^ f 2 

■npolyix)x dx 



(13) 



with Qb.o and Qm,o the cosmological density parameters of 
baryons and of total matter, respectively. 
asubsubsectionEntropy Injection Models 
X-ray observations of galaxy groups and clusters in- 
dicate that there may be an excess in gas entropy over 
that given by the polytropic model in the inner regions of 
low temperature cluster s and groups ( jPonman et al. I I1999I : 
iLlovd-Davies et al.|[200ol : IFinoguenov et al. 2007). This sug- 
gests that some non-gravitational processes, such as super- 
nova explosions and/or AGN feedback, may have raised 
the entropy of the gas. In our second model, we there- 
fore consider an entropy injection model similar to that in 
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iMoodlev etai] l|2008l ): IVoit et all l|2002l ). Throughout, we 
define the 'specific entropy' of the gas as 

S = T^n-^'\ (14) 

where rig = p^/diirip). Ignoring the details of the entropy 
injection process, we simply add a constant entropy term, 
5inj, to the polytropic entropy profile, Spoiy, which follows 
from Eq. (|14|l upon substitution of Eqs. ((6]) and 0. Defining 
-'^g(^) = ■^g(< r)/Mvir, with Afg the hot gas mass within r, 
we model the final specific entropy distribution of the gas as 



Sontl-Fg) — 5'poly(-Fg) + Sinj . 



(15) 



Based on X-ra y ob s ervations of clust ers (e.g. 
iLlovd-Davies et al.1 l2000l : IPonman etlo] Il999l ). we set 
the entropy floor to be either 100 keV cm^ (hereafter model 
'EntlOO') or 200keVcm^ (hereafter model 'Ent200'). The 
gas pressure is then Pcnt(r) oc S'cntPcnt , where pcnt is the 
gas density profile. In order to obtain pcnt{r), we solve the 
following equations, 



dPe, 



dr 
d£g 
dr 



— — Pcnt- 



GM{< r) 



— 4iYr 



2 Pent 



(16) 
(17) 



by setting the boundary conditions as follows. A natural 
boundary condition is given by -fg(O) = 0. However, the 
choice of the other boundary condition is less straightfor- 
ward. In the polytropic case described above, we have set 
fg(fvir) = (1 — /star) (Ohfi/^mfi)- Howcvcr, this couditiou 
is not expected to be valid in the presence of entropy in- 
jection, simply because the increase in entropy changes the 
temperature in the inner part of the halo, which may drive 
part of the gas out of the halo. Here, somewhat arbitrarily, 
we choose to set the boundary condition at r^ir by assuming 
that rcnt(''vir) = Tvir- At Icast this model has the desirable 
property that it reduces to the polytropic model described 



above in the limit Sin 



0. 



2.2.2 Physical Properties of the Hot Gas 

In order to link the above model to the SZ effect generated 
by hot electrons, we estimate the electron pressure as 



Po — ks ric Tc . 
The gas is assumed to be fully ionized, so that 
„ _ (1 + Xh) Pe 



(18) 



(19) 



with the hydrogen mass fraction xh = 0.76. The electron 
temperature is assumed to be the same as that of the gas, 
i.e., Te = Tg. 

Figd] shows the hot gas properties inside a halo of mass 
M = 10" /i~^M0 at redshift z = 0.1. The resuhs are 
shown for all three gas models discussed above: the poly- 
tropic model, EntlOO and Ent200. In all cases the mass 
fraction in stars and cold gas is assumed to be 0.1 times 
the total baryon fraction (i.e. /star = 0.1). Note that the 
polytropic model predicts a more concentrated gas profile 
and a higher electron pressure profile than the other two 
models. As we will see, this results in a stronger predicted 
total SZE signal. Entropy injection modifies the gas profile, 
resulting in a higher temperature in the inner regions of the 
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Figure 2. The ratio between the gas fraction in a halo and the 
universal gas fraction as a function of halo mass. Results are 
shown for the polytropic model (solid line), for which this ratio 
is equal to unity by construction, and for the entropy injection 
models EntlOO (dashed line) and Ent200 (dotted line). Stronger 
entropy injection results in a reduction of the hot gas density at 
small radii. 



halo. To satisfy the HE within the same dark matter halo, 
the gas density in the entropy injection models is therefore 
reduced relative to the polytropic case, and part of the gas 
now resides outside the (virial radius of the) halo. This is 
illustrated in Fig. [2l which shows the ratio between the gas 
fraction in a halo and the universal gas fraction, as a func- 
tion of halo mass for the three models considered here. In 
the case of a lO^** H'^Mq halo, this shows that entropy in- 
jection of lOOkeVcm"^ (200keVcm~'^) reduces the hot gas 
mass fraction inside the virial radius by ~ 15% (~ 30%), 
compared to the polytropic case without entropy injection. 

It is important to point out that the hot gas profile 
depends on the boundary conditions adopted. Current sim- 
ulations show that the gas tends to tr ace the dark mat- 
ter distribution in the outer part (e.g. Ijing fc Sutd l2000l : 
lYoshikawa et al.]|2000l : iLewis et aLll2000l ). but the tempera- 
ture profile of the hot gas is still poorly understood. Un- 
fortunately, the boundary conditions for the hot gas are 
not well constrained by observation, and different bound- 
ary conditions have been used in literature. For example, 
[Komatsu & Scljak (2001) derived a hot gas profile (KS 
model, hereafter) assuming polytropic gas and that the den- 
sity profile of the gas matches t hat of the dark matter in the 
outer regions of the halo, while lOstriker et al.l (|2005l ) set the 
boundary condition by requiring the gas surface pressure to 
match an exterior pressure. These different boundary con- 
ditions result in sizable differences in gas profile, especially 
in the outer parts of the halo. In our model, the hot gas 
distribution is assumed to have a sharp cutoff at about the 
virial radius. However, in reality, the gas in the vicinity of 
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Figure 1. Properties of tlie liot gas in a lialo with mass M = 10^* /i~^Mq at redsiiift z = 0.1. Results are sliown for all three gas 
models discussed in the text: the polytropic model, and the entropy injection models EntlOO and Ent200. In all three cases /star = 0.1. 
The various panels show the density (upper left-hand), temperature (upper right-hand), electron pressure (lower left-hand) and entropy 
(lower right-hand) of the gas, all as a function of halo-centric radius. 



dark matter halos may also be heated by the collapse of 
large scale structure, contributing to the SZE. Although we 
adopt the boundary conditions described above, we will oc- 
casionally comment on uncertainties that may arise from our 
specific choise of boundary conditions. 



3 THE EXPECTED SZE SIGNALS FROM 
DARK MATTER HALOS 

3.1 The SZ effect 

The temperature fluctuations caused by the SZE can be 
written as 



!MB 



(20) 
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where Tcmb ~ 2.73K l|Mather et al.l and x = 

hpv / fcBTcMB, with hp the Planck constant and v the pho- 
ton frequency. The first term on the right-hand side is the 
thermal SZE caused by the thermal motions of electrons; the 
second term is the kinematic SZE, caused by peculiar bulk 
motions along the line of sight of the object in question. The 
thermal SZE has a spectral shape given by 



+ 1 



(21) 



and a frequency-independent Compton parameter, j/comp, 
which is related to the projected electron pressure by: 



3 0"T 



(22) 



mcC^ 

with ar the Thomson cross section. In the isothermal case. 



with 



yc 



r = (Jt 



ricdl , 



(23) 



(24) 



the Thomson optical depth. The kinematic term is given by 
P — Vpcc/c, where Upoc is the bulk peculiar velocity along the 
line of sight. In what follows, we consider only the thermal 
part of the SZE, treating the kinematic part as 'contamina- 
tion'. 



3.2 The SZ effect produced by halos of a given 
mass 

Before exploring the SZ effect produced by a sample of 
groups that have similar masses, we first focus on the ther- 
mal SZE due to a single dark matter halo. Assuming spheri- 
cal symmetry, the Compton parameter j/comp around a halo 
of mass M can be written in terms of the 'projected electron 
pressure', Ep, as 



«/comp(-R|M) 



rricC 



(25) 



We can express Ep in terms of the cross-correlation between 
dark matter halos and the pressure profile of hot gas: 



Ep(i?|M) = 2ne / ^h.p 

JR 



(r|M) 



r dr 



R2 



(26) 



with R the projected distance to the halo center 

For convenience we will be working in Fourier space. 
The cross power spectrum between dark matter halos 
and the pressure field, Ph_p(fcjA'/), is related to the cross- 
correlation ^h.p(T'IAf) by 



Up{r\M) 



1 

2^ 



Ph,p(fciAf)^^fc'dfc. (27) 
kr 



This can be separated into two parts, a 1-halo term and a 
2-halo term. The 1-halo term describes the pressure profile 
of the halo's own hot gas atmosphere. For a halo of a given 
mass at a given redshift, this part of the power spectrum 
can be written as 



P^%(k\M) = ^^ip(fc|A'/), 

72c 



(28) 



where Up is the Fourier transform of the electron pressure 
profile, 

2 sin kr 



Up{k\M) = 



kr 



-n4r\M)Tcir\M)dr . (29) 



The 2-halo term describes the cross-correlation between halo 
centers and the pressure distribution of hot gas in other ha- 
los. We assume that on large scales the correlation function 
of halos is related to that of the dark matter by a linear bias 
relation. The 2-halo term can then be written as 

Ph,p(fc|M) = ^^^^b{M) n Up{k\M')n{M')b{M')dM' , 

(30) 

where Piin(fc) is the linear power spectr um of the density 
field, n{M) is the halo mass function fe.g. lPress &: Schechtej 
19741) and h(M) is the bias function for halos of mass M (e.g. 
Mo fc Wh ite 1996; Sheth et al. 2001). 

The SZE brightness is the change in the CMB specific 
intensity caused by the SZE, and can be written as 



A/(a;) = C,(x)Ioycomi, , 



where 



C(x) = 



2hp 



ksTc^ 



1 



(31) 



(32) 



(33) 



(e^ - 1)2 

and, as before, x = hpv / fcBlcMB. The integrated Compton 
parameter, Y , can be written as: 



Y : 



ycomp 

dA, 



(34) 



where dA — D^dfi, with Q, the solid angle of the halo and 
Da its comoving angular-diameter distance. In the isother- 
mal case, the value of Y is determined by the gas temper- 
ature and the total gas fraction, but is independent of the 
gas density profile. The SZ flux density, the change in the 
CMB flux density caused by the SZE within the solid angle 
S7, can then be written in terms of the integrated Compton 
parameter Y as 



2fep / fcBJcMB 



Y^x) 



(35) 



As one can see, the SZE brightness (flux) is proportional 
to the Compton j/comp {Y) parameter. The SZ brightness is 
redshift independent, but the total flux density of a cluster 
depends on its redshift through Da- 



4 CORRELATION BETWEEN GALAXIES 
AND HOT HALO GAS 

With a redshift survey of galaxies, another interesting quan- 
tity to study is the cross correlation between galaxies of 
different properties and hot halo gas. Since galaxy proper- 
ties are tightly correlated with the mass of their host dark 
matter halos, the measurement of such correlation can be 
used to constrain the hot gas distribution in dark matter 
halos as function of their mass. The correlation between the 
galaxy density and thermal S ZE has been studied in simula- 
tion and early WMAP data ljHernandez-Monteagudo et ahl 
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|2004 I2006I '). With upcoming SZE surveys, the correlation 
between galaxies and hot gas can be studied more de- 
tailedly. Here we develop a formalism to study the SZE 
around galaxies in a given luminosity bin. The interpre- 
tation of such observation does not require the identifi- 
cation of individual galaxy groups but is related to the 
distribution of galaxy luminosities as a function of halo 
mass. This connection between galaxy luminosity and halo 
mass is most conveniently described by the conditional 
luminosity function (hereafter CLF), ^{L\M)dL, which 
specifies the average number of galaxies with luminosities 
L ± d L/2 that reside in a halo of mass M l|Yang et al.l 
I2OO3I : Ivan den Bosch et all |2003| ). For a given cosmol- 
ogy, the CLF can be c onstrained us i ng galaxy clustering 
(e.g. lYang et all l2003l : ICooravl |2006|; Ivan den Bosch etlo] 
200j}_j_galaxy-galaxy lensing (e.g., iGuzik fc SeliakI l2002r 



Li et al. | ,2 009: Cacciato ct al. 200 9.) , sate llite kinematics 



fe.g.. Ivan den Bosch et al.ll2004lMore et a l. 2009, 201(f) and 
galaxy group catalogues (e.g., lYang et al.i ,2005, , .20081 ') . or 
any combination thereof. In this paper, we adopt the CLF 
parameterization motivated by the results obtained f rom a 
large galaxy group cat alogue (see Yang et "all |2008| ) with 
the parameters given in ICacciato et al.l ( 20091 ). which have 
been obtained using the combined constraints from the ob- 
served galaxy luminosity function, the luminosity depen- 
dence of galaxy clustering, and t he SPSS group catalog ue of 
Yang et al. (2007 ) . As shown in ICacciato et al.1 l|2009l ) and 



Li et al. ( 20091 ) the same CLF mo del also accurately matche s 



the galaxy-galaxy lensing data of lMandelbaum et al.l l|2006l ). 
Hence, in what follows we assume that the CLF is accurately 
known, and we focus on how to use it in order to extract in- 
formation regarding the hot gas properties in dark matter 
halos from the SZE-galaxy cross correlation. 

The expected SZE signal around a given galaxy also 
depends on the position of the galaxy within its host halo. 
For example, the SZE around central galaxies, which reside 
at the centers of their host halos, is expected to be quite 
different from that around satellite galaxies located at off- 
center positions. Hence, it is convenient to split the cross- 
power spectrum between galaxies and the hot gas pressure 
into 4 parts. The 1-halo central term, the 1-halo satellite 
term, the 2-halo central term, and the 2-halo satellite term: 



(36) 



P.Ak) = /c [Pg',^,'^(fc) + Pfp-(fe)] 

+ /s[j'g,'rW+^g,p"«]. 

where fc and fs are the central and satellite fractions, re- 
spectively, among among all galaxies in the luminosity bin 
being considered. Note that the power spectrum depends 
on both the redshift and luminosity of galaxies. For brevity, 
however, we will not write down this redshift dependence 
explicitly. 

In order to compute cross power spectrum of galaxies in 
a luminosity bin Li dL/2 we proceed as follows. Denote the 
probability that a central (satellite) galaxy with luminosity 
L resides in a halo of mass M by ■Px(M|L) [where x refers to 
either 'c' (central) or 's' (satellite)]. Using Bayes' theorem, 
we can write 



$x(L|M)n(M) 
ML) 



(37) 



where $x(i|M) is the CLF for central (satellite) galaxies in 



halos of mass M, and (f'xiL) is the corresponding luminosity 
function, which is related to the CLF according to 

ML)= ^4L\M)n{M) dM . (38) 
Jo 

Assuming that central galaxies always reside at the cen- 
ters of their host halos, we can write the 1-halo central term 



Pg!i"ik\L) = — Vc{M\L)up{k\M) dM . (39) 



Inserting Ea. (|37p into the above equation gives 
1 



P,i,^/(fc|L) 



nc4>c{L) Jo 



<l>c(-L|Af) Up{k\M) n{M) dM . 



(40) 

To calculate the 1-halo satellite term, one needs to know 
how satellites are distributed in their host halos. Here we 
make the assumption that satellite galaxies follow a number 
density distribution, Ua{r\M), that is similar to that of the 
dark matter particles; i.e., Us(r|M) oc (r/rs)~"'^(l -|-r/rs)~'^. 
The corresponding Fourier transform is 

Us(k\M) = 47r r™ uJr\M)^^^^^r''dr . (41) 
' kr 







The 1-halo satellite term can then be written as 



nc(t>s{L) Jo 



^4L\M)us{k\M)up{k\A.'I) n{M) dM . 

(42) 



The 2-halo term describes the correlation between 
galaxies and the hot gas pressure in halos other than their 
own host halos, and can be written as 



Pg'.^'"(fc|L) = Pnn{k)I4L)lM , 



(43) 



where Piin(fc) is the linear power-spectrum of the density 
field, and 



ML) 



<t>c{L) 



b{M) n{M) dM , (44) 



u4k\M)b{M)n{M)dM , (45) 



and 



Up{k\M)b{M)n{M)dM . 



(46) 



In practice, we consider the signal produced by galaxies 
in a finite luminosity bin, [Li, L2], which can be obtained by 
integrating the L-dependent quantities over L, and replacing 
$x(i|M) in the above equations by 



(iV)x(M) 



{L\M)dL. 



(47) 



5 DATA ISSUES 

In this section we discuss various characteristics of the on- 
going SZE surveys that are needed in order to make realistic 
predictions of the signal to be expected from the two dif- 
ferent analyses proposed here: stacking galaxy groups and 
cross-correlating galaxies with the SZ signal. 
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Table 2. Characteristics of Pla nck, SPT and ACT surveys. Dat a 
adapted from [Bartlctt ( 2006) H o et al.l (120091 ') and lSwetzl l l2009t) 



Figure 3. The redshift distribution for group of different halo 
masses in the group catalog. In the left panel, the solid line 
shows the redshift distribution of halos in the range [1012, 10^3] 
/i^^Mrri, and the dashed line shows that for halos with masses 



'i©, 
13 



ll3.5l 



h Mq. In the right panel. 



in [10", 10 

and dashed line represent halos of [lO^'^'^ , lO^*] 
[10", IQlS] h-lMg, respectively. 



the solid 



line 
and 



5.1 Group Catalogues and SZ Surveys 

One of the analyses we are proposing is to stack large num- 
bers of galaxy groups in order to probe the hot gas in rel- 
atively low mass haloes. Galaxy group catalogues are best 
obtained from large galaxy redshift surveys, and in what fol- 
lows we will focus on the SDSS. The wide sky coverage of 
the SDSS allows one to identify thousands of galaxy groups 
within a redshift of about 0.15 (see below). The typical angu- 
lar size of a 10^^ h~^yiQ halo is about 10 arcmin at z = 0.1. 
Hence, many of the SDSS groups are expected to have suf- 
ficiently large angular sizes to be resolved in future SZE 
surveys with typical resolutions of sub-arcmin to a few ar- 
cminutes. 

Various groups have c onstructed galaxy group cata- 



logues using t he SPSS fe.g. iGoto 



Berlind et al] l2006l : iKoester etal 



:ted galaxy 
l2005iriMiller 



et al 



20071 : lYang et al 



2005 



2005 



20071 ) . The most suita ble one for o ur pu rpose is the one 
published recently by lYang et al.l l|2007l . hereafter Y07). 
This group catalog (hereafter G CY07) was obtained from 
the S PSS Data Release 4 (PR4; lAdelman-McCarthv et all 
l200d)^ using th e halo- based adaptive group finder developed 
bv lYang et al.l (|2005l ). Each group is assigned a halo mass 
according to the total stellar mass it contains. As demon- 
strated in Y07, with the help of large mock redshift surveys, 
this group finder works well not only for rich groups but also 
for poor systems, facilitating studies of galaxy groups that 
cover a wide range in halo masses. 

In what follows, the GCY07 is used to demonstrate the 
feasibility of the analysis we are proposing here. In partic- 
ular, we use the GCY07 to estimate the number of groups 
in each mass bin and the corresponding observational noise 
level. 

FigO shows the redshift distribution of groups in differ- 
ent mass ranges. It is encouraging that the GCY07 provides 
a fairly large number of groups in the low-redshift Universe, 
which allows us to analyze groups in relatively narrow mass 
bins. For instance, there are 114 groups with masses between 
10"-^ and 10"-i h'^MQ in the redshift range 0.08 - 0.11. 
Table[T]lists the number of groups in different halo- mass and 



Name 


Freq. 


Res. FWHM 


Inst, noise 


Survey Area 




[GHz] 


[arcmin] 


[/iK/bcam] 


[dcg2] 


Planck 


143 


7.1 


6 


40000 




217 


5 


13 






353 


5 


40 




SPT 


150 


1 


10 


4000 




220 


0.7 


60 






275 


0.6 


100 




ACT 


148 


1.4 


15 


2000 




218 


1.3 








227 


0.9 







redshift bins. In the following sections, we will use the num- 
ber of groups in these samples as our basis for predicting 
the detectability of the SZE in upcoming surveys. 

Table [5] shows some important characteristics of the 
SPT, ACT and Planck surveys. For Planck, the coverage 
will be all-sky so that the number of groups available for 
stacking is limited by the optical survey. In this case, the 
numbers listed in Table [1] can be directly used for com- 
puting the predicted SZ signal. For the SPT survey, how- 
ever, the aimed coverage is about ~ 4000 square degrees. 
More importantly, it has zero overlap with the SPSS DR4. 
However, the 2dFGRS, which has a depth similar to the 
SDSS, has an overlap of about 1000 square degrees with the 
SPT. On-going optical surveys, such as the Dark Energy 
Survey (DES) Q and the next generation galaxy surveys, 
such as PanSTARR^ and LSSTQ will provide much larger 
and deeper optical samples covering all, or at least part of, 
the SPT survey area. In principle, the number of groups for 
stacking is then limited by the SPT survey area which is 
similar to that of the SDSS DR4. We would then expect the 
number of groups available in the same mass and redshift 
ranges to be similar to those listed in Table [T] However, 
these surveys will only yield photometric redshifts, which 
are far less reliable than the spectroscopic redshifts avail- 
able in, for example, the SDSS. Hence, group finders specif- 
ically designed for redshift s urveys , such as the halo-based 
group finder of lYang et all l|2005h used here, are not ex- 
pected to perform accurately on these surveys. However, in 
recent years some group finders have been developed speci- 
ficly for photometric surveys, and it has been demonstrated 
that they can achieve high completeness and purity (e.g. 
iKoester et"aIll2007l : [Milkeraitis et al.ll2010l ). Although these 
methods are mainly restricted to massive groups and clus- 
ters, this is not necessarily an important restriction for SZE 
studies such as those proposed here, which, as we demon- 
strate below, are only able to probe hot gas in relatively 
massive halos. In this work, we therefore make the optimistic 
assumption that future surveys will ultimately provide data 
in the SPT area of sufficient quality, so that our modeling 
based on the GCY07 is relevant. 



* https:/ /www. darkenergysurvey.org/ 
^ http:/ /pan-starrs. ifa.hawaii.edu/ 
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http:/ /www. lsst.org/ 
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Table 1. The number of halos (groups) in different mass bins and redshift bins in the group catalogue of Y07, which is based on the 
SDSS DR4. The first row lists the mass bins in log[A//( /i~-'Mq)]. The first column indicates the various redshift bins. 





[12.9,13.1] 


[13.4, 13.6] 


[13.9,14.1] 


[14.4, 14.6] 


[14.9, 15.1] 


[0.05,0.08] 


762 


252 


72 


18 





[0.08,0.11] 


1368 


448 


114 


16 





[0.11,0, 14] 


2836 


824 


169 


16 





[0.14,0.17] 


3707 


1208 


270 


38 


2 


[0.17,0.2] 


3864 


1594 


431 


77 


7 



For comparison, we also explore the detectability with 
the ACT survey. This survey has a beam size and instru- 
mental noise similar to the SPT, but a slightly smaller sur- 
vey area. If the corresponding optical survey has a similar 
depth, the expected number of groups in the ACT area will 
be about 50% of that in the SPT area. Thus, the signal to 
noise ratio expected from the ACT is about 70% of that 
from SPT. 



5.2 Noise and Contamination 

In order to examine the level of SZE that can be observed 
with the samples described above, we need to take account 
of noise and signal contamination in the observations. For 
the problem we are considering here, the main source of 
noise will be instrument noise, while there are three types of 
astrophysical sources that may cause signal contamination: 
(i) the primary CMB anisotropy, (ii) point sources, and (ill) 
the SZE produced by unresolved background clusters. We 
now discuss each of these in turn. 

The instrument noise (in /iK/beam) expected for the 
various SZ surveys is listed in Table [21 Throughout we will 
assume that the instrument noise for different galaxy groups 
is uncorrelated, so that stacking lowers the noise by a factor 
of y/N, with N the number of groups in the stack. 

The contamination by the primary CMB anisotropy 
is expected to dominate at large scales. The noise level 
due to the primary is about 100 fj.K per beam, which is 
much larger than the instrument noise. Fortunately, the 
primary anisotropy is frequency independent, while ther- 
mal SZE varies with frequency and vanishes at around 217 
GHz. Thus, the primary contamination can in principle be 
separated from the SZE by using multi-band observations. 
iPlagge et al] (|2010l ) studied the SZE profile of galaxy clus- 
ters in the SPT survey. They used the 220 GHz map to 
subtract the background contamination and produced a set 
of band-subtracted map with a depth less than 20 fiK. This 
provides hope that the contamination due to the primary 
CMB can be properly subtracted. In this paper, this con- 
tamination is ignored. 

Bright point sources such as quasars and star forming 
galaxies can contaminate the SZE map on small scales in 
the form of bright spots. However, since such sources are 
expected to be masked out, we do not consider them. Un- 
resolved point sources, on the other hand, can produce a sig- 
nificant contamination. In the case of unresolved IR point- 
sources the contamination leve l is expected to be compa - 
rable to the instrument noise () White fc Maiumdarlliooi l. 
Since the correlation among these sources is not expected 



to play a significant role (e.g. IWhite fc Maiumdarir2004l '). we 
can treat this contamination as un-correlated noise. 

Another class of unresolved point-sources that may 
be an important source of cont amination are radio galax- 
ies. Many investigation s (e.g. White fc Maiumdaj 120041 : 
IStaniszewski et ahl [2OO9I : IPlagge et al.l I2OI0I ) have shown 
that the contamination by unresolved radio sources is lower 
than that of the IR sources and unlikely to be an important 
contribution to the SZE noise. However, these sources may 
be correlated with the clusters and groups under investiga- 
tion. In principle, such correlation can be understood with 
observat ion in wav e -band s th at are not sensitiv e to S ZE. 
Recentlv lHall et al. I (!2010') and'Staniszewsk i et al.l (|2009l ) ar- 
gued that the clustering amplitude of radio galaxies is only 
a few percent of the mean background on arcminute scales. 
Thus, the contribution from clustered radio sources is at a 
level of a few tenth nK. Such noise is not important for 
groups with masses ~ 10^"* /i~^Mq or larger. However, since 
this noise does not decrease with stacking more groups, it 
may significantly affect results for groups of ~ 10^'^ h~^MQ 
and needs to be included. 

In this paper, we treat the point source contamination 
as un-correlated noise and account for it by simply doubling 
the instrument noise. 

Previous investigations have demonstrated that the 
contamination by the SZE background can be significant. 
For example, using light cones con structed from a n adi - 
abatic hydrodynamical simulation, iHallman et al.l ()2007l ) 
found that unresolved halos and filaments can contribute 
about 30% to the total flux of a SZE map, and that there 
is a significant chance (60%) for a single beam to contain 
multiple sources in a survey with a be am size as large a s 
that of the Planck survey. More recently, [Sh aw et al.l (|2008l ) 
studied this effect using a similar method. They derived a 
fitting formula for the SZE background fluctuations and es- 
timated their impact on the Y-M relation. They found that 
the SZE background contamination is about 10 to 30 percent 
of the total SZE fiux of a cluster at low redshift. In order 
to estimate the uncertainties introduced by this contamina- 
tion we compute the SZE angular power spectrum using our 
fiducial model (see Appendix), from which we construct a 
mock all-sky map of S ZE temperature flu ctuations, using the 
HEALPIX packagfQ l|G6rski et al.ll2005h . The angular reso- 
lution of our mock SZE map is 0.73arcmin^ per pixel (i.e., 
we set the HEALPIX parameter A'sido = 4096). Assuming 
that the mean contribution of the background sources can 
be subtracted, we can then estimate the rms of the fluctu- 
ations in the background on a given angular scale. On 1- 



http: / /healpix. jpl.nasa.gov 
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Figure 4. The Compton parameter j/comp as a function of the 
projected halo-centric distance. The contributions to the y pa- 
rameter are divided into 1-haIo term and 2-halo term. The three 
sets of Hues from high to low represent the results for halos at 
redshift z = 0.1 and with masses 10^^ h'^MQ, lO^* H'^Mq, and 
lO-'^ h~^MQ, respectively. 



0.01 I . . . . /••.-■.'.-I 

^q1,3 .qM ^q1 

M/(h"^Mo) 

Figure 5. The contribution of the 2-halo term to the total inte- 
grated Compton parameter, Y, within a projected radius equal 
to the halo virial radius as a function of halo mass for halos at 
z = 0.1. The solid, dashed and dotted lines are results for the 
fiducial model, EntlOO and Ent200, respectively. 



arcmin scale, we find an rms of 1.03 x 10"" arc min'^ for the 
Y par ameter, in agreement with the results of IShaw et al.l 
l|2008t ). Unlike the instrument noise, the noise generated by 
this contamination is expected to be correlated on small 
scale. To take such correlation into account, we mimic ob- 
servations separately for the SPT, ACT and Planck surveys. 
For Planck, we choose random directions in the mock sky 
and estimate the noise within solid angles chosen to match 
the virial radii of halos of different masses. For SPT and 
ACT, we estimate the noise in annuli around different ran- 
dom directions, with the sizes of the annuh chosen to match 
the radial bin sizes to be used. The SZE background fluctua- 
tion around clusters at two different directions is assumed to 
be uncorrelated. Hence, when stacking the signal from mul- 
tiple directions (i,e., multiple groups), the noise contribution 
due to this contamination source decreases with \/iV, where 
A'^ is the stacking number. This noise is added into our error 
budget by assuming that it is independent of the other noise 
sources. 

To summarize, our noise model consists of instrument 
noise, which we have artificially doubled to mimic the con- 
tribution of unresolved IR point sources and Radio point 
sources (assumed to be uncorrelated), plus the contamina- 
tion due to SZE background fluctuations obtained from our 
mock all-sky maps as described above. 



6 RESULTS 

6.1 Fiducial Model 

We first show the SZE for our fiducial model, where the gas 
is assumed to be in HE in NFW dark matter halos and to 
have a polytropic equation of state with F = 1.2. The cosmo- 
logical parameters a dopted are those from the WMAP3 data 
jSpergel et al.llioOTi '). Fig|4] shows the Compton parameter 
2/comp as a function of the projected halo-centric distance R. 
Results are shown for halos with three masses: 10", 10^"*, 
and 10^^ h~^MQ, respectively. The 1-halo and 2-halo terms 
are plotted separately using different line styles. While both 
terms decrease monotonically, the 1-halo term clearly dom- 
inates the SZE within the virial radius. The 2-halo terms 
for all the three halo masses have the same shape, and the 
difference in the amplitude is due to the mass dependence of 
the linear halo bias, b{M). For halos with M = 10"/i-^Mq, 
the two halo term becomes noticeable at the outer part of 
the halo, because the 1-halo term is relatively low. Note, 
however, that we have assumed that the hot halo gas only 
extends out to the virial radius. If the hot halo gas extends 
out to larger radii, the 1-halo term may still be important 
at larger radius. In particular, the models considered here 
ignore possible contribution from the warm-hot intergalactic 
medium (WHIM) associated with the filamentary and sheet- 
like structures in which the dark matter halos are embedded. 
As a simple test of the potential impact of such a WHIM 
component, we modified our fiducial model such that the 
hot gas profile extends to two times the virial radius. This 
boosts the SZE on scales around the virial radius by a fac- 
tor ~ 2, making it easily detectable for groups more massive 
than 10^'' h~^MQ. This suggests that the analysis proposed 
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here is also very promising in probing the relatively elusive 
WHIM in the direct vicinity of dark matter halos 

In FigOwe plot the fractional contribution of the 2-halo 
term to the total integrated Compton parameter Y within 
the projected radius R — r^ir as a function of halo mass. 
Results are shown for all three gas equations of state con- 
sidered here, as indicated. While the 2-halo contribution is 
negligible for the most massive groups (i.e., <, 5% for groups 
with M >, W^'^h"^ AIq), its contribution can reach as much 
as ~ 40% in haloes of ~ IO^/i'^Mq. As shown in Figdl en- 
tropy injection decreases the gas pressure inside dark matter 
halos, causing the amplitude of the SZE to decrease for both 
the 1-halo and 2-halo terms. Since the impact of entropy in- 
jection is more pronounced in less massive halos, and since 
the 2-halo term reflects the contribution from halos of all 
masses, entropy injection has a weaker effect on the 2-halo 
term than on the 1-halo term for low mass halos, while the 
opposite applies to massive halos. 



6.2 Dependence on Model Parameters 

Since many of the parameters in our model are still poorly 
constrained, we now investigate the effects on the predicted 
SZE signal from changing the following model parameters: 
the equation of state of the hot halo gas, the mass fraction 
of the hot gas, the concentration of dark matter halos, and 
cosmological parameters. For brevity, we only present re- 
sults for a halo with M = W^'^H'^Mq at redshift z = 0.1. 
FigE] shows the Compton parameter i/comp as a function of 
the projected halo-centric distance R. The top left panel 
shows the dependence on the gas equation of state. Here we 
compare the polytropic model with the two entropy injec- 
tion models, EntlOO and Ent200 (see Section[22]for details). 
Compared to the polytropic case, the gas pressure profiles of 
the entropy models are lower. This is easy to understand, as 
energy injection heats the gas at the halo center, reducing 
the density of the gas in the central regions. As expected, 
the reduction is more pronounced for a larger value of Sinj. 
Note that the 1-halo term in the entropy models extends to 
a slightly larger radius than in the polytropic model, which 
is due to the fact that the energy injection drives part of the 
hot gas out of the virial radius. The difference between the 
polytropic model and Ent200 is within a factor of two, and 
the effect is more important for lower-mass halos where the 
gravitational potential well is shallower. Note that the im- 
pact of an entropy floor is also evident at large radii, where 
the 2-halo term dominates, which is due to the additional 
heating reducing the hot gas pressure in individual halos. 

The top right panel of Fig[6] shows the model depen- 
dence on the hot gas fraction of the halo. As expected, j/comp 
is directly proportional to the gas fraction, making it trivial 
to scale our results for other choices of the gas fraction. 

Another model ingredient that remains somewhat un- 
certain are the halo concentration s. Observationally, these 
are only poorly constr ained (e.g. IComerford fc NataraianI 
120071 : lOguri et al.ll2009l '). which is why one typically resorts 
to the results from numerical simulations. These reveal a 
scatter in concentratio n paramet er, c, of about 30% at a 
given halo mass (e.g. Ijind |2000| ). More importantly, the 
mean halo concentra tion is found to dec r ease with increas- 
ing halo mass (e g. iBuUock et all I2001I: IZhao et al.ll2003l : 
iDolag et al.1l2004l : lMacci6 et al.ll2007l : IZhao et al.ll2009f ). but 



Table 3. Cosmological parameters used in the paper. 
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0-8 


WMAPl 
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0.7 


1.0 


0.9 


WMAP3 


0.238 


0.762 


0.041 


0.734 


0.951 
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the exact slope and normalization of this mass-dependence 
remain fairly uncertain, with different authors claiming re- 
lations that are significantly different. The middle left panel 
shows the effect of changing the concentrations of the dark 
matter halos by 20% (dotted lines) and 50% (dashed lines) 
relative to the fiducial value. Clearly, halo concentrations 
can have a significant impact on the SZE at small radii (and 
in the 2-halo term). Halos with larger (smaller) concentra- 
tions are more (less) centrally concentrated, which causes 
an increase (decrease) in the density of the hot gas at small 
radii. 

The middle right and bottom left panels of Fig. |B] show 
the effects of changing the cosmological parameters Qm and 
erg , respectively. Recall that for our fiducial model we adopt 
the WMAP3 cosmology (see TablegJ. In ACDM cosmologies, 
the 2-halo term, which reflects the clustering amplitude of 
halos, depends on both Qrn and ag through the halo mass 
function and halo bias function, while the 1-halo term de- 
pends on these parameters through the halo concentrations. 
Fig. [6] shows the impact on j/comp(J?) of 20% changes in the 
values of Qrn and erg relative to their fiducial values. As one 
can see, the SZE increases with as, while the dependence on 
Sim is rather weak. On scales dominated by the 1-halo term 
(i.e., R <, Tvir), the SZE depends only weakly on flm and erg. 
However, on larger scales, where the 2-halo term dominates, 
the effects are larger. In particular, an increase (decrease) 
of (jg by 20% results in an increase (decrease) of j/comp on 
large scales by a factor of ~ 2. A similar change in fim only 
affects i/comp at the 13% level. 



6.3 Predictions for the Planck Survey 

With the spatial resolution of Planck, low-mass halos are 
not spatially resolved, and ev en the highest mass h alos will 
only be marginally resolved (jAghanim eiaD[l993). In this 
case, observation can only be used to estimate the integrated 
Compton parameter Y of an entire group (or stack thereof). 
Fig[7] shows y as a function of group mass, M. Here Y is 
obtained by integrating j/comp within the virial radius of the 
halo. As discussed in Section O for a group with a given 
mass, M, the K-parameter decreases with the redshift of 
the group. For each mass bin we therefore average Y using 
the mass and redshift distributions of the groups in GCY07. 

To predict the error in a given mass bin, we use the 
number of groups in the GCY07 (see table [l| under the 
assumption that the noise for each group is independent of 
that of other groups. In that case the noise of the stacked 
signal decreases as N~^^^ , where A'^ is the number of groups 
in the stack. 

In the left-hand panel of Fig[7l the Y parameter is 
shown for our 3 different models for the gas equation of 
state. In all three cases, the Y parameter is almost the same 
at the high-mass end. For low mass groups, however, the 
entropy-injection models predict lower values for Y. For ex- 
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Figure 6. The model dependence of the Compton parameter j/comp- Here the halo has a mass of 10^ ?i~^M0 and is located at redshift 
z = 0.1. The top left panel shows the effect of changing the gas equation of state. The top right shows the effect of changing gas fraction. 
The middle left panel shows the effects of changing the halo concentration by 20% and 50% . The middle right and bottom left panels 
show the impact of changing the matter density parameter, Hm, and the power-spectrum normalization, erg, respectively, by 20% relative 
to the values for our fiducial (WMAP3) cosmology. 



ample, for halos with M ~ lO'^'^ h~^yiQ the Y value pre- 
dicted by the polytropic model is larger than that of the 
Ent200 model by a factor of about two. In the right-hand 
panel of FigO we show the effect of gas fraction. Here again 
the results are shown for gas fraction that is 2 and 0.5 times 
that adopted in the fiducial model. The larger error bars 
show the 3 <j detection limit for individual halos expected 
for Planck. The smaller error bars represent the uncertain- 



ties in the corresponding stacks of groups. Clearly, by stack- 
ing groups of similar masses, the total SZE flux around 
groups with masses M >^ 10^^'^ /i~^M0 can be constrained 
well. In particular, the GCY07 is large enough to constrain 
the hot gas fractions to an accuracy of <, 40% for groups 
with masses higher than 10^^'^ /I'^M©. However the sen- 
sitivity and spatial resolution of Planck is not sufficient to 
distinguish among the different gas equations of state con- 
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Figure 7. The integrated Compton parameter y as a function of iialo mass. Tiie halo redshifts are set to 2 = 0.1. The left panel shows 
the results assuming different gas equations of state, while the right panel shows the results assuming different gas fraction. The green 
errorbars represent the 3 cr instrument noise of the Planck telescope at 143 GHz, and the red errorbars shows the 3 cr error expected 
from the stack of groups in the corresponding mass range. 



sidered here, at least not when using a group catalogue the 
size of the GCY07. Larger and deeper optical surveys will 
be required to reduce the noise to sufficiently low levels. For 
example, if the number of groups is increased by a factor of 
4, the observation will be able to distinguish Ent200 and the 
fiducial model at ~ 3a confidence level. 



6.4 Predictions for the SPT and ACT Surveys 

In the case of the SPT, the resolution at 150 GHz is 1 arcmin 
with an instrument sensitivity of 10 /iK/beam. For compar- 
ison, the angular sizes corresponding to the virial radii of 
groups with masses of 10^^, 10^* and 10^^ h'^MQ at z = 0.1 
are 5.24, 11.3 and 24.3 arcmin, respectively. Thus, SPT has 
sufficient spatial resolution to probe the actual hot gas pro- 
files of SDSS groups. However, in order to obtain sufficient 
S/N, one needs to stack groups of similar masses together, 
especially at the low mass end. As for Planck, we base our 
analysis on galaxy groups in the GCY07, even though the 
area on the sky covered by this group catalogue (obtained 
from the SDSS DR4) has no overlap with the area covered 
by the SPT survey (see discussion in Section 15. ip . For a 
group of a given mass at a given redshift, we calculate the 
SZE profile and smooth it over a one-arcmin scale using a 
top-hat window. We then average the signal in each mass 
bin according to the host halo mass distribution given by 
the SDSS group catalog. 

The estimate is made in 11 logarithmic bins of halo- 
centric distance (R), from the halo center to a maximum 
of 10/i~^Mpc. To maximize the S/N, we stack groups in a 
wide redshift range, from 0.05 to 0.2, following the redshift 
distribution of the GCY07 group catalog. For each _R-bin, 



we compute the average signal as 



(48) 



where j/i,comp is the SZE signal of the ith group, and Wi is 
a weighting function chosen to be 



Wi = 



(49) 



with Di the angular diameter distance of the ith group. This 
weighting scheme is chosen to account for the fact that for 
haloes of the same mass, the SZE flux within an annulus of 
a certain angular distance is inversely proportional to the 
square of the distance of the halo. The noise of j/comp in a 
bin of R is estimated through 



ayiR) = 



T.,w,gI^^{R) 



(50) 



where ay^i{R) is the measurement noise of the ith group. 
Assuming that the noise contributions from different beams 
(i.e., different groups) are independent, we can estimate 
ay^i(R) through 



ay,i{R) 



(51) 



where (Tm is the measurement noise of t/comp per beam, and 
Ni{R) is the number of beams the annulus contains. Note 
that Ni{R) depends on i because the beam size corresponds 
to different real space area at different redshifts. The SZE 
background fluctuation, (Tbg, is estimated as described in 
Section 15.21 and is added in quadrature to ay to get the 
total uncertainty 



^tot = ^<(i?)+a2(i?). 



(52) 
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Figure 8. The average Compton parameter ycomp as a function of the projected halo-centric distance R. The top left panel shows the 
results of the fiducial model for different halo mass. The errorbars show the expected la uncertainty of the SPT for the corresponding 
stacks. The expected 1 a uncertainties of observing one single group is shown using blue solid line. The top right panel shows the results 
assuming different gas fraction. The bottom right panel shows the results assuming different gas equations of state, and the bottom left 
panel shows the result for different halo concentrations. In these three panels, we adopt a halo mass of 10" h-^MQ. The results are 
compared with SPT la detection limit (green error bars) for a single observation, and the detection limit for stacks (red error bars). We 
also plot the la detection limit of ACT(black error bars) for comparison. 



The upper left-hand panel of Fig[8] shows the aver- 
age Compton parameter, ^comp, as a function of the pro- 
jected halo-centric distance R for our fiducial model. Re- 
sults are shown for the three mass bins around 10^"', 10^'', 
and 10^^ /i~^Mq, respectively (see Table [!}. The expected 
1(7 uncertainties for the stacks of groups are shown with 



the errorbars. For comparison, the l-o uncertainty expected 
from a single group is indicated by the solid blue line. This 
shows that the SPT can only map the SZE profile of individ- 
ual halos with masses significantly larger than 10^'' /i~^Mq. 
However, by stacking the signal from groups in a catalogue 
the size of GCY07, one can study the (average) hot gas dis- 
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tribution around halos with masses as low as ~ 1O"^~^M0. 
The remaining three panels of Fig[8] show the predicted, 
average SZ profile for the stack of GCY07 groups in the 
10^* h~^MQ mass bin for different gas mass fractions (up- 
per right-hand panel), different gas equations of state (lower 
right-hand panel), and different halo concentrations (lower 
left-hand panel). Comparing these model predictions with 
the expected SPT sensitivity, it is clear that SPT can pro- 
vide stringent constraints on the properties of hot gas in 
dark matter halos. Even the 2-halo term can be detected 
with a group catalogue the size of GCY07, allowing one to 
probe hot gas in the infall regions around group-sized haloes 
(i.e., possibly associated with the WHIM inside filaments 
and pancakes). 

We also predict the signal and detection limit expected 
from ACT. The procedure of the calculation is the same as 
for SPT, except that the predicted signal is now smoothed 
with a 1.4 arcmin top-hat window to match the resolution of 
ACT at 148 GHz. Since the signals are binned into relative 
large annuli, the results expected from ACT and SPT are 
almost identical, except that the signal-to-noise is slightly 
lower for ACT. For comparison, we plot the error-bars ex- 
pected from the ACT in Fig|8] As one can see, ACT is still 
able to provide constraints on the hot gas fractions in halos 
with M > 10" /i-^Mq. 



6.4-1 The Impact of Mass Uncertainties 

In the analysis presented above, the calculations are made 
under the assumption that there are no errors in the halo 
mass assigned to each individual group. However, using 
mock galaxy redshift surveys, Y07 have shown that the error 
on the halo mass assigned to an individual group is of the 
order of 0.2 to 0.3 dex. To test the importantance of these 
errors, we repeat the same exersize as above, but this time 
adding a random 'error' to the halo mass of each group in 
the stack, drawn from a log-normal distribution with a rms 
of 0.2 dex. The 'perturbed' masses are then used to calcu- 
lated the SZE signal. Fig |12l shows a comparison between the 
results thus obtained with the 'perturbed' masses and those 
obtained assuming the halo masses are perfectly estimated. 
The difference caused by the mass uncertainties is negligi- 
ble at large radii (2-halo term). On small scales, where the 
1-halo term dominates, the mass errors cause an increase 
in j/comp of about 7%, and the 1-halo term now extends to 
slightly larger radii. This is due to the inclusion of massive 
halos in the tail of the mass distribution. 

In real observation, this effect may be quantified by ap- 
plying the same analysis to mock group catalogs selected in 
the same way as the real catalog. Using these mock cata- 
logs, one can estimate the uncertainties in the assigned halo 
masses, which in turn can be incorporated in the modeling, 
similar to what we have done here. 



6.5 The Cross Correlation between Galaxies and 
the Compton Parameter 

Another way to probe the properties of the hot gas in dark 
matter halos is to study how the hot gas correlates with 
galaxies of different luminosity. One advantage of such anal- 
ysis is that the number of galaxies that can be used is 




Figure 12. The effect of mass assignment errors. The solid line in 
the left panel shows the original mass distribution in the mass bin 
[10"-9, 10" l] /i-iM0 as given by GCY07, while the dashed line 
shows the distribution including the broadening modeled with a 
lognormal distribution with a dispersion of 0.2 dex. The right 
panel shows the results obtained with these two mass distribu- 
tions. 



large, so that one can stack the SZE signals around many 
galaxies, thus achieving high signal-to-noise. Even more im- 
portantly, measurement of the galaxy-SZE cross correlation 
does not require a selection of galaxy groups, which always 
carries some uncertainties and which typically requires spec- 
troscopic redshifts for individual galaxies. In fact, as we 
demonstrate below, the precision of photometric redshift 
may be sufficient to obtain reliable measurements of the 
galaxy-SZE cross correlation. Since photometric redshifts 
are much easier to obtain than their spectroscopic coun- 
terparts, the galaxy-SZE correlation can be used to study 
the hot gas evolution out to higher redshift. The stacking 
results can then be interpreted with the CLE discussed in 
Section [3] The solid lines in the upper two panels of Fig. [5] 
show ycomp{R) around stacks of galaxies with absolute r- 
band magnitudes in the ranges —20 < Mr — 5\ogh ^ —19 
(left-hand panels) and -21.5 < Mr - blogh -22 (right- 
hand panels). The dotted, dashed, and dot-dashed curves 
show the contributions from the 1-halo central, the 1-halo 
satellite and the 2-halo terms, respectively. These are com- 
puted using t he method descr i bed in Section|3]with the CLE 
parameters of lCacciato et all (|2009l ) . Errorbars are obtained 
using the same method as described in Section 16.41 except 
that rather than using the numbers of galaxy groups in the 
GCY07, we use the numbers of galaxies in the spectroscopic 
sample of the SDSS DR4 with absolute r-band magnitudes 
in the respective bins, and with redshifts 0.05 ^ z ^ 0.20. 
The error-bars are estimated using the total independent 
angular area of all the angular annuli (corresponding to the 
bin size in R) around all galaxies in the corresponding lu- 
minosity bin, in comparison to the SPT (ACT) beam size. 
The lower panels of Fig. |5]show the host halo mass distribu- 
tions for central galaxies (solid lines) and satellite galaxies 
(dashed lines) in the corresponding magnitude bins. For clar- 
ity, the distributions are individually normalized, while the 
satellite fractions in each luminosity bin are indicated in the 
corresponding panels. 

As one can see, the SZE signal for the brighter sample 
is dominated by different contributions on different scales. 
The 1-halo central term dominates the signal up to a scale of 
about 0.3/i"^Mpc. On scales of 0.3/i"^Mpc <R< 2/i"^Mpc 
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Figure 9. The predicted SZE around galaxies in two different luminosity bins (as indicated by the r-band absolute magnitudes in the 
low two panels). In the upper panels, the predicted j/comp parameter is plotted as a function of the projected distance to the galaxy. 
The SZE is decoinijosed into a 1-lialo central term (dotted lines), a. 1-lialo satellite term (dashed lines), and a 2- halo term (dash-dotted 
lines). The total signal is shown by the solid line. The Itr expected uncertainties from SPT are shown as the errorbars on solid circles. 
We also plot the la confidence level expected from the ACT survey as error-bars on triangles. In the two lower panels, the host halo 
mass distribution of the galaxies used to obtain the results shown in the upper two panels are shown for central galaxies (solid line) and 
satellites (dashed line). The distributions are normalized individually. The satellite fraction /gat in each luminosity bin is marked in the 
respective panel. 



the signal is dominated by the 1-halo sateUite term, while 
the 2-halo term dominates on larger scales. For the fainter 

sample, the situation is difTcrcnt. Since the host halos 
of these central galaxies are mostly low mass halos with 
10";i-^M© < M < 10^^ h-^MQ (see lower left-hand panel), 



the 1-halo-central term extends only to relatively small 
scales and never dominates the signal. The 1-halo satel- 
lite term dominates the SZ signal from the center out to 
~ 1 /i~^Mpc, after which the 2-halo term dominates. This 
is important, as it indicates that even very accurate mea- 
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Figure 10. The figure shows how the haloes in different mass range contribute to the SZE signal around galaxies that in certain luminosity 
bins. The absolute magnitude range of galaxies are marked on top right of each subplot. Different linestyles represent contribution from 
halo of different mass range as shown on top left panel. The unit is in h~^'MQ 



surements of the average 3/comp(^?) around relatively faint 
galaxies yields virtually no constraints on the properties of 
hot gas in low mass haloes (i.e., the haloes that host central 
galaxies of those luminosities). Rather, since the SZ signal 
scales strongly with halo mass {Y oc Af° with a ~ 1.6), the 
signal on small scales is dominated by the 1-halo satellites 
term, even though satellites only make up a relatively small 
fraction (~ 21%) of all galaxies in that luminosity bin. 

This is also evident from Fig. [I0]which shows how galax- 
ies in halos of different masses contribute to the SZE effect. 
Results are shown for four different luminosity bins, as in- 



dicated. In the brightest bin, [-22.5,-21.5], halos in the mass 
range 10^^ - 10^^ h~^MQ dominate the signal in the inner 
part, while more massive halos only dominate the signal be- 
tween 0.4 and lh~^ Mpc. This can be understood with the 
help of Fig[51 From the inner to outer parts, the signal of 
the brightest bin is first dominated by the 1-halo central 
term and then by the 1-halo satellite term. Since host ha- 
los of satellites on average are more massive than those of 
central galaxies of the same luminosity, the 1-halo satellite 
term dominates on intermediate scales. For faint luminosity 
bins, the SZ signal on small scales ( ^ 1 /i~^Mpc) is always 
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Figure 11. The model dependence of j/comp(^?) predicted around galaxies in the luminosity bin [—22,-21.5]. The Ic expected uncer- 
tainties from SPT are shown as the errorbars on solid circles. Icr detection limit for ACT survey using error bar on triangle. The upper 
left panel shows the results for different equations of state. The upper right panel shows the results for different gas fractions. The lower 
left panel shows the results for two different set of cosmological parameters. And the lower right panel shows the effect of changing the 
halo concentration. 



dominated by the most massive halos. Although only a rela- 
tively small fraction of faint galaxies reside in these massive 
halos (as satellites) their contribution to the SZ effect is 
larger than that from the far more numerous centrals in less 
massive halos. 

Fig. I 111 shows the predictions for galaxies in the [-22.5,- 
21.5] luminosity bin for different gas equations of state (up- 
per left-hand panel), different gas mass fractions (upper 
right-hand panel), different cosmological parameters (lower 



left-hand panel), and different halo concentrations (lower 
right-hand panel). Here again we see that the statistical 
uncertainties expected from a survey like SPT are much 
smaller than the differences between different models, in- 
dicating that an analysis along these lines can put tight 
constraints on the hot gas properties in dark matter halos 
spanning a relatively wide range in masses. 

The lower left-hand panel of Fig. [11] shows a compari- 
son between ycomp(-R) expected around a stack of galaxies in 
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the [-22.5,-21.5] luminosity bin in the WMAPl and WMAP3 
cosmologies (see Table [3] for the corresponding cosmological 
parameters). In addition to changing the cosmological pa- 
rameters here we have also changed the CLF that is used to 
predict UcompiR)- For each cos mology, the best-fit CL F pa- 
rameters have been obtained bv lCacciato et al] (|2009t ) using 
the combined constraints from the observed galaxy luminos- 
ity function, the luminosity depend ence of galaxy clu ster- 



ing, and the SPSS group catal ogue of lYang et al.l (120071 ). As 
shown in lCacciato et all (|2009| ). both models fit the observed 
abundance and clustering of galax;ies equally well. How- 
ever, whereas the WMAP3 mo del simultaneously matche s 
the galaxy-galaxy lensing data of lMandelbaum et al.l (|2006h , 
the WMAPl clearly overpredicts the lensing signal (i.e., the 
WMAPl model predicts mass-to-light ratios that are too 
high). Hence, galaxy-galaxy lensing can be used to break the 
degeneracy b etwe en cosmolo gy and halo occ upation statis- 
tics (see also I Yang et al.ll2006l : iLi et al.ll2009l ). Interestingly, 
the results presented here suggest that the cross correlation 
between galaxies and the Compton parameter can also be 
used to discriminate between these different models. How- 
ever, it is clear from Fig. [11] that a change in cosmologi- 
cal parameters has a similar effect as a change in the hot 
gas fraction or a change in halo concentration. Hence, it 
will be important to have precise constraints on cosmologi- 
cal parameters in order to use the galax;y-Compton param- 
eter cross correlation to constrain the properties of hot halo 
gcis. The precision of upcoming CMB observations, such as 
Planck, will be able to pin down all important cosmological 
parameters to very high accuracy, making the SZE analysis 
presented here a powerful tool to study the gas distribution 
in dark matter halos. 

Finally, for comparison, we also plot the error-bars ex- 
pected from the ACT survey. Clearly, the uncertainties in 
ACT are slightly larger than that of SPT and can also tightly 
constrain the properties of the hot gas in halos that host 
these galaxies. 



6.5.1 Photometric Redshifts 

When calculating the cross-correlation between galaxies and 
hot gas, we have assumed that the redshifts of the galaxies 
have no errors. Although this is a reasonable assumption to 
make when using spectroscopic surveys, we now investigate 
how significantly redshift errors, such as those present in 
a photometric redshift survey, impact on the analysis. Red- 
shift errors introduce errors in both the absolute magnitudes 
(luminosities) and the projected distances. In order to esti- 
mate the resulting impact on ycomp(^Z) we add a redshift 
error (and its associated luminosity error) to each galaxy in 
the [-21.5, -22] luminosity bin of the SDSS DR4 catalog, 
using a Gaussian distribution with standard deviations of 
(l + z) X 5% and (1-l-z) x 3%. This roughly mimics the red- 
shift errors expected from future photometric surveys, such 
as the LSST. Next we compute the galaxy- Compton param- 
eter cross correlation using the new redshifts. Results are 
shown in Fig. 1131 The change in the cross correlation due to 
photo-z errors of 5% (3%) is negligible at small radii, but 
increases to ~ 20% (~ 10%) at 1 - 2 h'^Mpc. We also find 
that, if one can exclude the 20 percentile of the galaxies with 
the largest redshift errors, the measurement errors caused by 
the photo-z errors shrink to about half this value. This sug- 
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Figure 13. This figure shows how the photo-z error affects the 
observed j/comp- The solid Une shows the Compton parameter 
l/comp(J2) in the luminosity bin [-22,-21.5]. The dashed and dotted 
Unes show the results of the calculations that include a photo-z 
error of 5% and 3% in (H-z), respectively. 



gests that the cross correlation analysis considered here can 
be extended to photometric surveys, but that precise model- 
ing of the galaxy-Compton parameter cross correlation may 
require some treatment of the photo-z error distribution. 



6.6 Model based on observed X-ray gas profile 

In previous sections, the hot gas profile is modeled based 
on the assumption of hydrostatic equilibrium between the 
hot gas and the dark matter halo. This model is roughly 
consistent with the observe d X-ray luminoisty an d gas frac- 
tion of galaxy c lusters (e.g. Moodlev et al.ll200§ ). However, 
recent work bv lKomatsu et ah I (|2010l ) showed that the gas 
profile thus deri ved is different from that derived from X- 
ray observation ijArnaud et al.l[20ld . AlO, hereafter). Since 
our fiducial model is very similar to the KS model (see Sec- 
tion l2.2.2)) . there is the same uncertainty due to the assumed 
hot gas density profile. 

In Fig |14l we compare our hot gas pressure profile with 
that derived from AlO, with the la tter obtained u si ng th e 
formula given in the Appendix of iKomatsu et aP (|2010h . 
The horizontal axis is scaled by -R500 , the radii within which 
the mean dens i ty is 500 times the critical density. As in 
iKomatsu et al.l ()2010h . we also find that the pressure profile 
predicted by our model is flatter and more extended. This 
difference causes a significant difference in the predicted SZE 
at large radii. 

One possibility to alleviate the tension between the two 
models is to change the gas fraction and boundary condition 
in our model. In our model, the boundary condition is set by 
T{r vii) — Tvir, but X-ray observations indicate that the tem- 
perature drops more rapidly in the outer part of the halo. For 



20 Ran Li et. al 



0.01 001: 




0.1 1.0 

R/Rsoo 



Figure 14. The hot gas profile of o ur fiducial model (so lid lines) 
in comparison with that derived in l lArnaud et al.ll2010h (dashed 
lines). The upper panel shows the hot gas pressure profile as 
a function of radius. The bottom panel shows the Compton y- 
parameter profile. In each panel, the upper set of lines are for 
the halo of 10^* /i^^Mq and the lower set of lines for halo of 
10^3 H-'^Mq. The horizontal axis shows radii scaled by the cor- 
responding i?500- 



some halo massed, the slope of the X-ray derived pressure 
profile can be fitted by decreasing the temperature at virial 
radius. However, in hydrostatic equilibrium the low bound- 
ary temperature results in a higher inner gas density, and 
the amplitude of SZE at small scales becomes higher than 
that based on X-ray observations. For a halo of 10^'*/i~^Mq, 
we can decrease the discrepancy between our fiducial pro- 
file and X-ray profile by setting T(rvir) = O.STvir and the 
gas to total fraction /gas = 0.08. However, we cannot find a 
single boundary condition that can fit the X-ray profile for 
all halo masses. Because of these uncertainties, the bound- 
ary condition should be treated as a free parameter when 
studying the hot gas properties of groups with future SZE 
observation. 

Observationally, X-ray observations lack the sensitiv- 
ity to reliably probe hot gas beyond -R500, because the X- 
ray luminosity scales with the square of the gas density. In 
A09, the pressure pro file is derived using REXCESS data 
ijBohringer et al.ll2007l ). where X-ray measurements are lim- 
ited within J?5oo. Thus, the pressure profile at larger radii is 
still uncertain. 

In order to test how the uncertainties in the hot gas 
pressure profile affect our prediction of the SZE, we com- 
pare in Fig llSI the cross-correlation between the SZE and 
galaxies predicted by our fiducial model with that predicted 
by the X-ray observation-based model. The predictions of 
our fiducial model is a factor of 2 to 3 larger. This indicates 
that the observation of the cross-correlation between SZE 
and galaxies can be used to distinguish these two models. 




fiducial 

A09 

10'° I I 10'° I 

0.1 1.0 10.0 0.1 1.0 10.0 

R/(h''Mpc) R/(h''Mpc) 

Figure 15. The j/comp(R) around galaxies in luminosity bins 
[-22, -21.5] (left) and [-20, -19] (right) predicted by our fiducial 
model (solid lines) in comparison with th at predicted with the 
formula derived from llArnaud et al.|[201Cll ) (dashed lines). 



7 CONCLUSIONS 

In this paper we have used simplified models to demonstrate 
the potential of using the SZE to probe the hot gas expected 
to be associated with dark matter halos. We have shown 
that by stacking SZE maps expected from the Planck, ACT 
and/or SPT surveys, one can probe the hot gas properties in 
galaxy groups with halo masses down to ~ 10^^ h~^yiQ. The 
SZE for halos with similar masses are examined in terms of 
halo models. Splitting the total SZE signal into 1-halo and 
2-halo terms, we have shown that in high-mass halos the 
1-halo term dominates the signal within the virial radius; in 
low-mass halos the 2-halo term is also significant on small 
scales. The model predictions are sensitive to the amount of 
hot gas assumed to be in halos, and our results show that 
a SPT-like survey can provide a stringent constraint on the 
hot gas fraction for halos with masses M >, 10^^ /i~^M0. 
Furthermore, such observations can also provide stringent 
constraints on the equation of state of the hot gas as well as 

the concentrations of dark matter halos. 

Using WMAP-7-year data, iKomatsu et al.l (|2010l ) 
stacked clusters in X-ray catalog and compared the 
resulti ng SZE profile with t he prediction of the KS 
model (|Komatsu fc Selia^ 120011 ) that uses a polytropic gas 
equation of state very similar to the polytropic model con- 
sidered here. They found that this model over-predicts the 
SZE by a factor o f 1.5 t o 2. In a recent study using SPT 
data, iLueker et al.l (|2009l ) also found that the observed SZE 
power spectrum is lower than that expected from polytropic 
model. However, as we have shown above, entropy injection 
can reduce the amplitude of the SZE by a factor of 1.5 at 
the inner parts of halos with masses lO^*/i~^M0, suggesting 
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that non-gravitational heating may have played an impor- 
tant role. 

We have also explored the idea of using the cross corre- 
lation between hot gas and galaxies of different luminosity 
to probe the hot gas in dark matter halos with the use of 
the CLF. Since the number of galaxies that can be used 
is large, especially in the forthcoming large photometric 
surveys where accurate photometric redshifts of individual 
galax;ies can be obtained, the measurement can be made to 
high precision. Our results show that, with the help of CLF 
modeling, one can constrain the hot gas profile in halos with 
masses down to 10^^ h~^MQ. Cosmological parameters, es- 
pecially the value of ug, can affect galaxy - hot gas cross- 
correlation. Thus, in order to use the observed SZE to con- 
strain the hot gas properties in dark matter halos, precise 
constraints on cosmological parameters are required. 

In summary, the upcoming SZE surveys are expected to 
provide an important avenue to study the hot gas properties 
in dark matter halos. Combining with large optical galaxy 
surveys and using the stacking method proposed here, the 
SZE to be observed in forthcoming surveys will allow us to 
study in great detail the hot gas properties in galaxy groups 
where most galaxies reside. This, in turn, will shed new light 
on how galaxy form and evolve in dark matter halos. 
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APPENDIX A: ANGULAR POWER 
SPECTRUM AND SZE MAP 

We write the SZE power spectrum, Psz(fe, z), as 

z) = P^^ik, z) + (fc, z) , (Al) 

where the 1-halo contribution is 

P,'^^{k,z) = lt^{k,z), (A2) 

and the 2-halo contribution is 

P^^{k,z) = Pii4k,z)lls,. (A3) 

Here 

l!Mk,z)= n{M,z)ul{k\M,z)dM . (A4) 

and 

/2,sz(A;,z)= / n{M,z)up{k\M,z)b{M,z)dM . (A5) 

ObservationaJly, one typically measures the angular power 
spectrum, which is the integration of the power spectrum 
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along the line-of-sight. To derive the angular power spec- 
trum, we expand the CMB temperature fluctuations due to 
the SZE in spherical harmonics: 



AT 

Tc 



MB 



AT 

Tc 



MB 



The angular power spectrum can then be written as: 
AT \ / AT 

im \TCMB . 



CTSll'Sr, 



(A6) 



(A7) 



Using the Limber approximation ()Limberlll953h . we have 

where d?V / AzAQ, is the comoving volume per unit redshift 
per unit solid angle, and the projected power spectrum, 
Pi,sz{l,z), is related to the 3-dimensional power spectrum 
via 



[ll.Bzil, Z) + /2,sz(^, Z)] 



with 



(7T 



(A9) 



(AlO) 



Here Dc{z) is the comoving distance from redshift z, and 
a{z) is the scale factor. 

We carry out the integration using Eo. dASIl from red- 
shift Zmin = to Zmax ~ 5. For the mass limits, we use 
Mmin = 5 X 10" /i-^M© and M^ax = 5 x 10^^ /i-^Mg. 
We then use the angular power-spectrum as a n input to the 
synfa st subroutine in the HEALPIX package (|G6rski et alj 
I2OO5I ) in order to generate an all sky SZE map. In FiglAl 
we show the standard deviation of the background SZE flux 
measured within an angular radius 6 (solid line) and com- 
pare it with the instrumental noise levels of SPT, ACT and 
Planck. The background noise is much smaller than the in- 
strumental noise. We also compare o ur result with the fitting 
formula given by Sha w et al.l (|2008l ) (dotted line) who use 
the background SZE noise measured from a light cone sim- 
ulation with cosmological parameters similar to ours. Note 
the good agreement between the two. 




4 6 
angular raclius(6) 

Figure Al. Tlie fluctuation of tiie SZE background as a function 
of angular radius 9. The solid lin e shows our r esult, and the dotted 
line shows the fitting formula of lShaw et al. I ll2008l 'l. The triangles 
show the instrumental noise levels of SPT (150GHz), ACT trian- 
gles show the instrumental noise levels of SPT (143GHz), ACT 
(148 GHz) and Planck (150GHz), respectively. 
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